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An initial study of the extremely high energy (EHE) physics capability of the IceCube neutrino observatory 
is demonstrated by considering a GZK mechanism neutrino production model, which is a guaranteed source for 
EHE neutrinos. We study EHE event properties in the energy range 10 5 < E < 10 11 GeV observed by IceCube 
using detailed Monte Carlo simulation. Results of a simulation study show that about 0.7 EHE neutrino events 
per year is expected with the full IceCube configuration over a 0.03 atmospheric muon background which passed 
an EHE event criteria. It is also shown that with the present partial IceCube detector with the same criteria is 
capable of studying EHE physics and the sensitivity improves with the number of deployed strings in the ice. 



1. Introduction 

It is well known that there exists extremely 
high energy (EHE) particles in the universe with 
energies approaching up to ~ 10 20 eV [1]. These 
EHE cosmic rays (EHECRs) may produce neu- 
trinos by various mechanisms, namely when they 
interact with surrounding matter or photon fields. 
The neutrinos are generated by the decays of ir 
mesons vi&ir^ — > /x^z^t — > e^feV^ process. In the 
EHE region, collisions of EHECRs and the cosmic 
microwave background photons produce cosmo- 
genic neutrinos [2], a consequence of a process 
known as the Greisen-Zatsepin-Kuzmin (GZK) 
mechanism [3]. 

In general neutrinos are unique probes for ex- 
ploring the high energy universe mainly because 
of two distinct features. The first feature is that 
the existence of cosmic neutrinos implies that 
there are hadron beams because only high en- 
ergy hadronic processes can produce neutrinos. 
The second feature is that neutrinos can pene- 
trate dense matter and radiation fields because of 
their small cross-sections, and therefore can prop- 
agate over cosmological distances. This implies 
that by observing neutrinos, one is able to learn 
about stars and galaxies which are surrounded 
by dense gases or energetic objects at very large 
distances. Such information is obscured when ob- 
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serving with photons (e.g., radio, visible light, X- 
ray and 7-ray). Because of these features, fur- 
thermore, detection of EHE neutrinos may shed 
light onto one of the most difficult questions in 
modern astrophysics, how and where are EHE- 
CRs produced. 

It is generally very difficult to detect high en- 
ergy neutrinos because one needs a large target 
volume for neutrinos to interact and produce a 
detectable signal in a reasonable time scale. The 
IceCube neutrino observatory, however, provides 
a rare opportunity to overcome this difficulty. Its 
km 3 fiducial volume which uses clean glacier ice 
deep below the surface at the South Pole is a pow- 
erful tool to search for EHE neutrinos. In this 
paper we discuss the capability of the IceCube 
experiment to detect EHE neutrinos. 

The paper is outlined as follows: In the next 
section, the overview of the IceCube observatory 
and its present status are briefly mentioned. The 
section 3 describes the characteristics of EHE 
neutrino events in the IceCube detector. We ex- 
plain signals and backgrounds in our EHE neu- 
trino search and a possible strategy to discrimi- 
nate neutrino events by considering their event 
topologies. In the section 4, the main results 
from the detailed Monte Carlo simulation study 
are shown, referring to the preliminary numbers 
concerning effective area and event rate for cos- 
mogenic EHE neutrinos. We summarize our con- 
clusion and future prospects at the end. 
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2. The IceCube Neutrino Observatory 

IceCube is a next-generation cubic-kilometer 
scale high energy cosmic neutrino telescope cur- 
rently under construction and in operation at the 
geographic South Pole. It uses 3 km thick glacial 
ice as a Cherenkov medium. Cherenkov photons 
emitted from relativistic charged particles such 
as muons are received by an array of Digital Op- 
tical Modules (DOMs) which amplify and digi- 
tally sample the high-speed photomultiplicr tube 
(PMT) pulses in situ. Each DOM encloses 10" 
R7081-02 PMT made by Hamamatsu Photonics 
in a transparent pressure sphere along with the 
high voltage system, a LED flasher board for opti- 
cal calibration in ice, and a digital readout board. 
The deep-ice DOMs are deployed along electrical 
cable bundles which carry power and information 
between the DOMs and surface electronics. The 
cable assemblies, often called strings, are dropped 
into holes drilled to a depth of 2450 meters. The 
DOMs occupy the last 1000 meters at intervals 
of 17 meters where the glacial ice is transparent. 
DOMs are also frozen into tanks located at the 
surface near the top of each hole which consti- 
tutes an air shower array called IceTop. IceTop 
provides us with the capability to study the atmo- 
spheric muon background reliably. This feature 
plays a key role in the IceCube EHE neutrino 
search as we will mention later. The strings and 
tanks are arranged in a hexagonal lattice pattern 
with a spacing of approximately 125 meters. At 
completion (planned to be 2011), the array will 
comprise 4200 in-ice DOMs on 70 strings and 320 
modules in the surface array. Currently, IceCube 
includes 9 in-ice strings (540 DOMs) and 32 Ice- 
Top tanks (68 DOMs). 

Events recorded by IceCube are generally cate- 
gorized by geometrical patterns of photon distri- 
butions from the minimum ionizing in-ice charged 
particles. The "track" events are initiated ei- 
ther from penetrating neutrino- induced muons or 
taus, or from muons coming from extensive air 
shower cascades above the surface i.e., the at- 
mospheric muons. Directional reconstruction of 
these tracks suppresses the atmospheric muon 
background. By selecting neutrino events that 
are up going, this effectively eliminates the at- 



mospheric muon background since they are down 
going. The other event pattern known as a cas- 
cade" is induced by an electromagnetic (EM) cas- 
cade via the v e charged-current interaction or 
hadronic cascades via the neutral-current inter- 
action of all neutrino flavors. Those cascades 
generated inside IceCube detection volume emit 
spherical shower-like patterns of Cherenkov light. 
Because muon and tau tracks have a longer prop- 
agation length in matter/ice than that of the elec- 
tron and hadron cascades, tracks generally have 
a larger effective area. 

EHE events that are expected to be seen by the 
IceCube observatory can be considered beloinging 
to neither class, however, because of difference in 
their major energy-loss processes. At the high- 
est energies, photon emission patterns differ and 
their remarkably sizable energy deposit in the Ice- 
Cube array may create a different characteristic 
signature. We describe this EHE event topology 
in the next section. 

3. EHE neutrino events in IceCube 

In the EHE region, because of the increase 
of the neutrino cross-section with energy, neu- 
trinos are more likely to be involved in interac- 
tions with matter during their propagation than 
to penetrate through the Earth. Charged lep- 
tons and hadrons are generated in these inter- 
actions and the secondary produced /x's and r's 
travel the Earth losing their energies by under- 
going many radiative reactions, i.e. EM cascades 
generated by e ± pair creation, Bremsstrahlung, 
and hadronic cascades generated by the photonu- 
clear interactions. The IceCube detector is to ob- 
serve these secondary /z's and r's as a main de- 
tectable channel of EHE signals [4]. There are 
two prominent characteristics for these EHE neu- 
trino events. One is that a major fraction of them 
arrive at the detector with a down-going geome- 
try because the mean free path of EHE neutrinos 
is much shorter than the typical path length in 
the Earth. Another is that the muon/tau "track" 
is accompanied with many "cascades" originat- 
ing from the various radiative energy loss pro- 
cesses. The first characteristic makes it difficult 
to discriminate between the atmospheric muon 
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background events and neutrino induced events 
by using the Earth as a filter. The second point, 
however, provides another way to distinguish the 
EHE signals from background because the energy 
loss due to stochastic radiative processes is pro- 
portional to the energy of muon(tau). 

A measurement of the energy deposit inside 
the IceCube detection volume, therefore, leads 
to an estimation of the track energy. The es- 
timated track energy tells if it is of cosmic ori- 
gin. This is because the expected spectra of sec- 
ondary /i's and t's generated from the GZK neu- 
trinos is much harder than that of atmospheric 
muons [4], and the measured energy (or its indi- 
cator) should be able to exclude the atmospheric 
muon events in a relatively straightforward man- 
ner. Figure 1 shows examples of simulated Ice- 
Cube events at different energies. One can see 
that the EHE muon radiates a large number of 
Cherenkov photons. The amount of Cherenkov 
photons in an event, the "brightness" of the event 
in other words, appears to indicate the track en- 
ergy. Detailed Monte Carlo study has shown that 
the number of detected photons is a robust energy 
indicator to search for EHE signals, which is de- 
scribed in the following section. 
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Figure 1. Examples of simulated IceCube events. 
Left panel shows a 100 TeV muon track repre- 
senting the conventional IceCube event while the 
right panel indicates a 9 EeV muon EHE event. 
Colored circles denote DOMs with more than one 
photon, the size of the circles represents the num- 
ber of photons and the color indicate the first pho- 
ton arrival time at each DOMs. x, y and z axes 
are distances from the central position of the Ice- 
Cube detector. 



4. Results 

4.1. Simulation Setup and Signal Domain 
Criteria 

Transportation of the EHE neutrinos and their 
secondary particles from the Earth's surface to 
the detector depth is calculated by the JULIeT 
package [5] as described in Ref. [4]. The IceCube 
Monte Carlo simulation package then generates 
events with energies and intensities following the 
obtained fluxes. The primary cosmic neutrino 
flux in this paper is assumed to be GZK cos- 
mogenic neutrinos as calculated in Ref. [6]. For 
the atmospheric muon flux, which is considered 
our main background, we take the analytically 
fitted result using a Corsika simulation [7] assum- 
ing an initial cosmic-ray proton flux of the form of 
E~ 3 . We should note that the atmospheric muon 
background in the EHE regime is highly uncer- 
tain due to our poor knowledge of the EHECR 
mass composition and the muon bundle inten- 



sity. The background intensities will be deduced 
from the IceCube data in the future. The present 
simulation chain considers events induced by the 
charged leptons that interact within 860m from 
the center of detector because the secondary pro- 
duced muons and taus are our main EHE events 
in IceCube [4]. The energy range under consid- 
eration is between 10 5 and 10 11 GeV near the 
detector. 

As described in Section 3, the total amount of 
Cherenkov light detected for an event is expected 
to be sensitive to the energy of incoming EHE 
muon or tau. The IceCube detector can mea- 
sure the Cherenkov luminosity per event as the 
event-total number of photo-electrons (NPE) de- 
tected by each DOM. Figure 2 shows distributions 
of NPE as a function of the muon and tau en- 
ergy. It shows a correlation between NPE and the 
muon/tau energy going through the IceCube vol- 
ume. Notice that the track geometry also affects 
the observed NPEs since a track passing further 
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Figure 2. Simulated event distributions on a 
plane of NPE and the charged lepton energy. En- 
ergy plotted here is the value when the incoming 
particle is at the surface of IceCube detection vol- 
ume which is defined as a sphere of 860 meters in 
radius in the present study. The left plot shows 
muon events and the right plot shows tau events. 
A suppression of energy loss for taus compared 
to that of muons and the contributions from tau- 
decays are visible in the right plot. 



away from the DOMs will result in less Cherenkov 
light within the IceCube detector. A distinct fea- 
ture in the plots is that the primary energy of 
the incoming leptons is still a deciding factor in 
the observed NPEs, which implies that NPE is 
a robust indicator of the lepton energy without 
having to rely on event geometry reconstruction. 

The zenith angle distributions are quite differ- 
ent for muon and taus originating from the prop- 
agation of cosmogenic GZK neutrinos and atmo- 
spheric muons, however. As down-going GZK 
muons and taus have a higher tendency to ar- 
rive from near the horizontal direction than the 
atmospheric background [4] , one can expect a dif- 
ference in the distribution of zenith angles and 
NPE for EHE signals and background. Plotted in 
Figure 3 are the simulated event distributions on 
the plane of NPE and zenith angle of muons and 
taus expected from the GZK model and the atmo- 
spheric muon backgrounds, which behave consis- 
tently with our expectation. Taking advantage of 
such differences, we introduce the signal domain 
on the NPE-Zenith angle plane where the GZK 



Figure 3. Event distribution on the plane of NPE 
and cosine of zenith angle obtained by Monte 
Carlo simulation. Plotted in the left and mid- 
dle are that of GZK neutrino-induced muon and 
tau signals respectively and the background at- 
mospheric muons on the right. 



events dominates over the background defined as 
5 x 10 5 ifcos0>O.l, 



NPE > 



I 1 x 10° otherwise. 



The expected event rate and effective area with 
this criteria are discussed next. 

4.2. Event Rate 

The expected signal and background event rate 
are shown in Figure 4 as functions of NPE and 
zenith angle. 

Major parts of the dominant atmospheric back- 
ground in the left plots are rejected by the selec- 
tion criteria defined as Eq. (1) in the right plots. 
With these selection criterion, the GZK muon and 
tau signal events dominate. It is also found that 
the number of tau events are comparable to or 
even greater than that of muon events. In the 
EHE regime, the tau does not decay but is sub- 
ject to radiative energy loss processes. Its heav- 
ier mass suppresses the energy loss compared to 
muons with the same energy, which gives a higher 
survival probability from near the horizontal di- 
rection. 

Shown in Table 1 are the event rates from dif- 
ferent string numbers deployed in ice. The same 
selection criteria are applied for all the string con- 
figurations, including the present 9 string setup. 
Event rate increases almost constantly as more 
strings are installed. 
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Figure 4. /i/r GZK signal and atmospheric /x 
event rate as function of NPE (upper) and cos (9) 
(lower). Blue and black markers indicate GZK 
/i, r respectively and red markers indicate atmo- 
spheric fi. Left plots display event rates without 
background cut and right plots show after the sig- 
nal domain selection. Multiple muon events are 
not included in the atmospheric muon flux. 



4.3. Effective Area 

The effective area resulted from the EHE sig- 
nal domain cut is plotted in Figure 5 as a func- 
tion of incoming charged lepton energy for dif- 
ferent zenith angles. The NPE-based EHE event 
selection criteria decreases effective area in lower 
energy region while it enhances the detection ef- 
ficiency for events with energies above 10 s GeV. 
These events are so luminous that the DOMs re- 
ceive many photons even from the tracks propa- 
gating far outside the instrumented detector vol- 
ume. The area reaches ~ 2 km 2 , a factor of two 
larger than the IceCube physical cross section in 
energy region above f0 10 GeV. The area for up- 
going geometry is slightly larger as the NPE cut 
is loose for those geometries. 



Table 1 

The preliminary IceCube EHE event rate for 
different deployed in-ice string configurations. 
Muons and taus produced from the GZK model 
and the muon atmospheric background are 
shown. 



string number 9 20 40 60 80 

GZK fj, 0.067 012 021 028 0.35 

GZK r 0.055 O.f 1 0.f9 0.25 0.31 

atmospheric ^ 0.009 0.010 0.012 0.025 0.033 
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Figure 5. Muon (left) and tau (right) effec- 
tive area as a function of incoming track energy 
with the 80 string configuration based on the 
EHE signal criteria. Dashed line denotes horizon- 
tal events and dotted line represents down-going 
events. Solid line indicates up-going events, al- 
though the probability of having EHE neutrinos 
with up-going geometry is very low, because of 
the Earth's sheltering effect. 



5. Summary and Future Prospective 

We have considered the capability of the Ice- 
Cube detector in the search for EHE neutrinos 
without extensive energy reconstruction. The 
first level background rejection and signal selec- 
tion is performed in terms of NPE and MC truth 
incoming direction. With this selection criteria, it 
is shown by a Monte Carlo simulation that about 
0.7 EHE GZK events per year is expected with 
the full IceCube configuration over a background 
of 0.03 atmospheric muons. 

The expected GZK event rate with the present 
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9 string configuration is about 0.12 events per 
year with the same event selection. The num- 
ber increases with the number of strings deployed 
from 9 strings to the full number of strings. 

The study presented in this paper has not 
involved energy and geometry reconstruction of 
muon and tau tracks. We expect that good ge- 
ometry reconstructions will support our results 
based on the Monte Carlo truth angle informa- 
tion. In addition, utilizing energy reconstruc- 
tion should lead to a significant improvement of 
detector performance, especially for uncontaincd 
events that traverse outside the IceCube instru- 
mentation volume. Ongoing development of EHE 
reconstruction methods will improve our sensitiv- 
ity to EHE cosmic neutrinos in near future. 

The results we have shown here heavily rely 
on NPE-Energy relations. These contain sizable 
uncertainties due to our incomplete understand- 
ing of the ice properties and detector responses 
to high NPE signals. The ice properties include 
the optical properties of the glacier ice and the 
behavior of the Cherenkov photon propagation in 
the hole ice which was melted and refrozen during 
the detector deployment. In order to reduce the 
systematic errors in the present study, the abso- 
lutely calibrated light source called the standard 
candle (SC) has been deployed in neighborhood 
of the DOMs in ice. The SC consists of a nitro- 
gen laser, the calibrated PMT, and the related 
optics system and is able to emit pulses which 
are approximately equivalent to a PeV EM cas- 
cade. Analyzing the data generated by the SC 
should reduce current uncertainties and improve 
the reliability of the EHE neutrino search. 

The intensity of the atmospheric muon back- 
ground in the EHE energy region is not well un- 
derstood as we mentioned in Section 4, and we 
have assumed the Corsika-based analytical model 
in the present simulation study. In the future, 
however, it is expected that reasonable statistics 
of the data with NPE below ~ 10 5 , allows us 
to build a reasonable empirical model to predict 
number of the background events in the EHE sig- 
nal domain. Wc also would like to point out that 
in any case our sensitivity would not be changed 
significantly by varying the definition of energy 
based background rejection factor because the 



GZK neutrino induced muons and taus exhibit 
an extremely hard energy distribution in contrast 
to the very soft (~ E~ 4 ) spectrum of the atmo- 
spheric muons. Besides, the atmospheric muon 
events can be experimentally identified and ex- 
cluded by looking for coincidence events between 
the deep-ice IceCube strings and the IceTop sur- 
face array. With this power of the IceTop air 
shower array, the hybrid analysis using the in- 
ice DOMs and the surface DOMs, even with a 
limited data statistics, provides us further infor- 
mation about the background, including events 
consist of multiple muons, in the very high en- 
ergy region. 
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